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Abstract 

We discuss three different scenarios recently proposed to account for the non-Fermi hquid behavior near antifer- 
romagnetic (AFM) quantum critical points in heavy-Fermion systems: (i) scattering of Fermi liquid quasiparticles 
by strong spin fluctuations near the spin-density-wave instability, (ii) the breakdown of the Kondo effect due to 
the competition with the RKKY interaction, and (iii) the formation of magnetic regions due to rare conflgurations 
of the disorder. Here we focus on the first scenario and show that it explains in some detail the anomalous tem- 
perature dependence of the resistivity observed e.g. in CePd2Si2, CeNi2Ge2 or Celua. The interplay of strongly 
anisotropic scattering due to critical spin-fluctuations and weak isotropic impurity scattering leads to a regime 
with a resistivity p ~ po + Ap(r), Ap{T) = T'^^'^f{T/po) cx for sufficiently large T and smaU po- 
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1. Introduction 

In recent years an increasing number of heavy- 
fermion metals near an AFM quantum critical 
point (QCP) were shown to display striking devia- 
tions from conventional Fermi-liquid behavior [1- 
3]. In particular, the resistivity rises as a function 
of temperature with exponents smaller than 2, the 
specific heat diverges logarithmically or shows a 
Vt cusp, and the susceptibility shows anomalous 
corrections of the form , with a < 1 . 

The common feature of the systems we are 
considering is that, as a function of some control 
parameter, which can be doping, pressure or mag- 
netic field, antiferromagnetic order is suppressed. 
The anomalous behavior mentioned above oc- 
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curs in the vicinity of the QCP defined by the 
vanishing of the Neel temperature [4]. We stress 
that both magnetic and non-magnetic phases are 
heavy-Fermion metals which display, for exam- 
ple, the characteristically large specific heat. It is 
commonly believed that the QCP results from the 
competition between the Kondo screening of mag- 
netic moments and the AFM correlations induced 
by the RKKY interaction [5]. At least three dif- 
ferent theoretical scenarios have been proposed to 
account for the observed behavior near the QCP. 

The first scenario is based on the the assumption 
that in a heavy Fermion system below a scale Tk 
(see Fig. 1) the low energy excitations are (heavy) 
quasi particles and their collective excitations. 
In this case the QCP should be in the same uni- 
versality class as the weak-coupling spin-density 
wave (SDW) transition in a Fermi liquid studied 
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Fig. 1. Scenario 1: Spin-density-wave transition in a Fermi 
liquid. In a Kondo lattice, below a scale T^, heavy quasi 

particles describe well the low energy excitations of the 
Fermi liquid. They undergo a AFM phase transition which 
is in the universality class of a SDW transition [6]. The 
QPT is driven e.g. by pressure p, doping or magnetic fields. 
Right figure: Near the QCP, scattering on the Fermi surface 
is enhanced along "hot lines" with ejj fa ek±Q « /i, where 
Q is the ordering vector of the AFM. In the remaining 
"cold regions" inelastic scattering is weak. 
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Fig. 2. Scenario 2: Breakdown of the Fermi liquid. The 
Kondo effect, which tries to screen local moments, breaks 
down near the QCP due to the competition with the RKKY 
interaction and singular magnetic fluctuations [9,11]. An 
effective Kondo temperature vanishes at the QCP. 

by Hertz [6]. More precisely, the behavior near 
the QCP is determined by the mutual interaction 
between Landau damped spin fluctuations and 
inclastically scattered quasi-particles. The bulk of 
this paper concerns the transport properties near 
the QCP, an issue of great experimental relevance. 
In spite of the substantial amount of work in 
this area [7] the interplay of strongly anisotropic 
inelastic scattering and weak isotropic impurity 
scattering had not been properly analyzed until 
recently [8]. 

While this first scenario appears internally con- 
sistent, it is far from clear that the notion of heavy 
quasi particles survives near the QCP. For exam- 
ple, one can envision a situation where the Kondo 



Fig. 3. Scenario 3: Disorder. Due to disorder, the effec- 
tive Neel temperature varies within a sample. Near the 
QCP these fluctuations get more and more important and 
change the critical behavior. In a region around the QCP, 
the thermodynamics is dominated by rare disorder config- 
urations which lead to the formation of magnetic domains 
even in the paramagnetic phase. 



effect breaks down near the QCP as a result of the 
competition of the Kondo effect with the RKKY 
interaction and singular magnetic fluctuations. In 
this second scenario the effective "Kondo" temper- 
ature vanishes at the QCP as shown in Fig. 2 
and the quantum critical behavior is dominated by 
local magnetic fluctuations. A microscopic justifi- 
cation of this scenario is still missing but might be 
realized along the lines proposed recently by Q. Si 
et al [9]. At present, CeCue-xAux seems to be the 
best experimental candidate for such a scenario 
[10,11]. At the QCP, the dynamic susceptibility has 
been fitted to x(q,w)-i « f {q) + a{bT - ico)'^ [10]. 
Astonishingly, the anomalous exponent a w 0.8 
shows up not only near the ordering wave vector 
Q of the AFM, but over the entire Brillouin zone. 
This suggests a local origin for the non- Fermi liq- 
uid behavior as might be induced by a breakdown 
of the Kondo effect, consistent with the ideas men- 
tioned above. 

The third scenario I would like to mention is 
driven by the presence of sufficiently strong disor- 
der. Of course, if typical fluctuations of the effective 
Neel temperature due to disorder are large com- 
pared to the distance of the quantum critical point 
(Harris critcrium), disorder effects arc important. 
Perhaps more importantly, even at some distance 
form the QCP, in the non-magnetic phase, a rare 
configuration of impurities can lead to a small mag- 
netically ordered region. Such exponentially rare 
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effects may dominate some of the thermodynamic 

properties of the system in a finite region around 
the QCP [14] (Fig. 3); e.g. cy/T oc x oc r-i+^, 
where the exponent A < 1 should vary continu- 
ously as a function of the distance of the QCP. Sev- 
eral strongly disordered uranium compounds can 
be fitted to this scenario [14] at least over some 
range. An interesting systematic approach to this 
problem was suggested by Narayanan et al [13], 
who studied the transition in a double e expansion 
in 4 + e space and et time dimensions. They find 
that rare events are indeed relevant in the rcnor- 
malization group (RG) sense if the number of or- 
der parameter components is smaller than 4. The 
flow to strong disorder in the RG equations can be 
interpreted as the signature of the onset of a Grif- 
fiths phase, but also other strong-coupling fixed 
points might be possible. It would be interesting to 
investigate whether the Kondo eflfect, which is not 
captured in an expansion around time dimension 
will qualitatively change the RG fiow. 



2. Resistivity 

In the following we will concentrate on the trans- 
port properties within the first scenario. We con- 
sider a situation where the influence of disorder on 
the magnetism is small, but we will keep track of 
the effects of weak impiuity scattering. 

We propose [8] to study the interplay of isotropic 
elastic scattering and anisotropic scattering from 
spin waves within a semi-classical Boltzmann equa- 
tion treatment of electrons interacting with spin- 
fluctuations and impurities. We use the fact that 
the former are described by a theory above the up- 
per critical dimension [6] and, due to the ohmic 
damping of the magnetic excitations in a metal, is 
characterized by a dynamical exponent ^; = 2. As 
a result, the spin-fluctuation spectrum in the para- 
magnetic phase can be modeled by [6,7,15,16] 



where qo and T are characteristic momentum and 

energy scales; Q is the ordering wave vector in the 
AFM phase and £, is the AFM correlation length 
which is near the QCP of the form [6] 
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r measures the distance from the quantum critical 
point and is e.g. proportional top — pc in an exper- 
iment tuned by pressure p. For the purposes of our 
numerical calculations we set c = 1, even at the 
QCP for r = the temperature dependence of £, 
does not influence the low temperature resistivity. 

Our starting point is the Boltzmann equation 
with a quasi-particle distribution function /k = 

— <i?k(9/k/9ek) linearized around the Fermi dis- 
tribution with a collision term 



coll 1^/ 



r 



(3) 

Imxk-k'(ek-ek') 



Here ^^^p and are transition rates for isotropic 
impurity scattering and inelastic scattering from 
spin fluctuations, respectively, and n° is the Bose 
function. In the following we measure the resistiv- 
ity in imits of pM = 3?i,g|/(7re^?;|,) and express 
our results in terms of the dimensionless elastic 
scattering rate x = irgf^^/ {2g'g) [18], the reduced 
temperature t = T/T and the effective distance 
from the critical point r (Eq. 1). In these units the 
residual resistivity po is given by po = xpM- Eq- 
(3) tacitly assumes that the spin-fluctuations stay 
in equilibrium, an approximation valid if the spin 
fluctuations can loose their momentum effectively 
by Umklapp or impurity scattering. 

At T = the scattering is purely elastic and 
the quasiparticle distribution function is given by 



$k(T = 0) = eEvk/(<7LpiVF). For T > 0, 



Xq(w) 



i/W + (q±Q)V9o-i'^/r' 



(1) 



$0 

the AFM spin fluctuations are strongest near the 
Q vectors in reciprocal space where the AFM or- 
der develops [16]. Accordingly, the electrons with 
an energy ek ~ A* will scatter strongly from the 
spin fluctuations near "hot lines" on the Fermi sur- 
face where Ck ~ Ck-FQ ~ M (right part of Fig. 1). 
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For sufficiently small temperatures, the distribu- 
tion function will change only in the neighborhood 
of the hot lines as long as the scattering rates in 
the "cold regions" (Fig. 1) are dominated by im- 
purities, i.e. for t"^ <C x. The strong scattering due 
to spin-fluctuations will equilibrate the points k 
and k ± Q along the hot lines on the Fermi sur- 
face and the non-equilibrium distribution function 
is written in the form <I>k ~ (1 — Pit)^'^ +-Pk*I'k±Q' 
where pk approaches 1/2 near the hot lines and 
vanishes further away. It is helpful to use a coor- 
dinate system, where vectors on the Fermi surface 
k = k// -I- kj^ are split up in a vector kjj on the 
hot line and a perpendicular vector kj^ _L k//, Vk. 
After rescaling the momenta using k = KqoVi the 
Boltzmann equation in the limit t < \fx can be 
rewritten in the following way 

t j dn' {pk +PkH±ci+K' - ^)I{rlt + M^k') 

+CkHa;pK = (4) 



I{y) = J dzz^n%z){l + n\z))/{y^ + z% 



(5) 



where k' = k',, + k'j_, Mm^i 
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and Ck„ = {2iTfNFVi,„+Q/ {-Kql). From (4) it is 
clear that only the combinations t/x and t/r deter- 
mine Ph. Accordingly, the resistivity obeys a scal- 
ing relation for t <C ^/x <C 1, r <C 1: 



Ap d t 2 r 2 

«t2/ , 

Pm \ X X 



(6) 
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,r <t <x^ 

r 2 _ 5-d _ 

, max[x'*-i , \/xr * \ <t < y/x 



/r 2 J f < niin[r, y/xr " ] 



where for completeness we have given the result in 
d dimension. Note, however, that the Boltzmann 
equation approach presumably breaks down md = 
2 as the "hot" electrons acquire a lifetime a VT. 
An overview of the behavior of Ap in = 3 is given 
in Fig. 4. 



nonuniversal regime 
t > V5. p^t^ 




Fig. 4. In the scaling limit t,x,r — > 0, t/x,r/x — > const. 
the resistivity is "universal" Ap/pM = t^'''^ f{t/x, r /x) for 
t < \/x, r < 1, where t oc T, x oc po and r oc p — pc 
measure the temperature, the amount of disorder and the 
distance from the QCP. The plot shows the qualitative be- 
havior of the scaling function / in the various regimes. For 
min[a;, ^/rx < t < ^/x\ the resistivity rises linearly with 
temperature. Thermodynamic quantities show a crossover 
to Fermi liquid behavior at the scale t = r, while in trans- 
port a behavior is recovered at a much smaller scale 
miii[a;, \/rx]. The dashed line serves as a reminder that at 
lowest tciupcratuics effects which are not included in our 
approach become important, e.g. interference effects of dis- 
order and interactions [12] or a disorder induced change of 
the spin fluctuation spectrum [13]. 

The asymptotic behavior can be easily extracted 
from (4) by realizing, that pk varies smoothly 
parallel to the hot lines, allowing us to replace 
Pkjj±Q-FK' = Pk'^^k'^ byp^ii^'^ in (4) and to perform 
the integration. For dominant impurity scatter- 
ing, t < min[x, \/rx], one can use a perturbative 
expansion in the first term in (4) and obtains 



Pm 



E 



dk, 



((Vk|| - Vk||±Qjn)2 



h - 



(vk|| X Vk||±Qje||| 

t 

OC 

I J yn\y){l + n°{y))Im^j+iy 



(7) 



( 7rC(3/2)r(5/2) 



,t>r 



(8) 



2V2 

/^V^r ,t<rT 
where en is a unit vector parallel to the hot line and 
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Fig. 5. Scaling plot of /(t/x, 0) = Ap{t/x) /t^^'^ at the 
QCP (r = 0) using a numerical solution of the Boltzmann 
equations for various amounts of disorder x. The straight 
solid lines show the limiting cases Ap oc t^/^ for t <ti x and 
Ap oc t for X <^ t <^ y/x. Deviations from scaling are large 
for t/x > 1/y/x. 

sums over all hot lines. Only in this dirty limit, 
the commonly used [15] approximation is valid to 
calculate the resistivity from the average of the 
scattering rate. 

In the limit ^/x > t > min[a;, ^/rx\, the resistiv- 
ity rises linearly with temperature 



Ap 

Pm 



dkii 



((vk|| - Vk||±Q,)n)' 
|(vk|| X Vk||±Qje|| 



, 5/4 

,3/4 Vk||±Q,/vk|| 



where the only slightly varying function 5° — 
J^^dnp^icLja) w TT is calculated from the solu- 
tion of the two coupled one-dimensional integral 
equations 



Pk 



dn'- 



dn' 



(Pk+Pk' - l)(l-a^) 
(n'^/a + K'^a - 2KK'a)3/2 

fe+PK' - 1)(1-Q^) 
(n'^a + K''^/a - 2KK'a)3/2 

±1 are solved by — 







(9) 



which in the limit a 

l/{{l+a)+K^/2y^). 

For 1 ^ t 3> y/x the resistivity is much less 
universal and as pointed out by Hlubina and Rice 
[7] p/pm ~ i^, where the prefactor depends on 
all details of the scattering processes in the "cold 
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Fig. 6. Resistivity as a function of T^'^ in CePd2Si2 taken 
from [2] (left figure) compared to our calculation (right fig- 
ure) for X = 0.01. The insets show the corresponding loga- 
rithmic derivative of p(T)-p(0). Below ft! 400mK CePd2Si2 
is superconducting (lower inset). Note the offset of the line 
T = in both plots. 

regions" . The full scaling curve for r = and the 
deviations from scaling are shown Fig. 5. 



3. Conclusions 

A main qualitative difference of the SDW sce- 
nario (Fig. 1) discussed above and a scenario, 
where the Kondo effect breaks down at the transi- 
tion (Fig. 2), is that in the first case only a small 
fraction of the quasi particles near the "hot lines" 
is strongly scattered by the critical fluctuations 
(Fig. 1) while in the second case the full Fermi 
surface is affected. 

At first glance, there is a simple way to distin- 
guish these two situations [11]: for an ultra-clean 
Fermi liquid interacting with spin-fluctuations the 
quasi particles with the longest lifetime short- 
circuit those with a strong scattering rate and de- 
termine the resistivity. Accordingly, one expects a 
Fermi liquid like resistivity in a system without 
disorder within the SDW scenario [7] (the limit 
X <^t^ discussed above). This would suggest, 
that a low-temperature resistivity of the form T" 
with a < 2 in a clean system signalizes a total 
breakdown of the Fermi liquid. Our calculations, 
however, show that the ultra-clean limit is irrele- 
vant even for samples with a RRR of 10'^ or 10"' [8] 
because very weak impurity scattering leads to the 
qualitatively different resistivity discussed above. 
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Our results imply that in all situations where two 
scattering mechanisms compete, one anisotropic 
and affecting only parts of the Fermi surface, the 
other isotropic, the transport is quite subtle and 
the interplay of the two mechanisms can lead 
to qualitatively new effects. For example, Math- 
iessen's rule is not valid in this regime. 

Taking this into account, it seems that quite a 
number of heavy Fermion systems are good candi- 
dates to be described within our first scenario, a 
SDW transition in a Fermi Hquid (e.g. CePd2Si2, 
CeNisGes, Celng, CeCusSiz, CeNiGaa [2,3]). The 
qualitative trends are well described by the theory: 
Ap cx T^'^ in the more dirty systems and Ap oc T" 
with a close to 1 in the cleaner systems. Also the 
cleaner systems seem to show anomalous behav- 
ior in the resistivity (but not in thermodynamics) 
even at some distance away from the QCP, con- 
sistent with our results (see Fig. 4). In Fig. 6 the 
resistivity of a quite clean sample of CePd2Si2 [2] 
is compared to a typical solution of the Boltzmann 
equation for x — 0.01 (see [8] for the precise model 
used). It is important to point out, that not only a 
similar effective exponent shows up in both theory 
and experiment, but that it is also observed over a 
similar range O.lpo < ^p{T) < lOpo- Certainly, it 
has to be checked carefully whether and to what ex- 
tend other predictions of the SDW scenario are ob- 
served experimentally. E.g. the predicted [6] pres- 
sure dependence of the Neel temperature Tjv cx; 
(Pc—p)'^^'^ appears to disagree with the experimen- 
tally observed linear behavior in CePd2Si2 [2]. 

The best test of our prediction would be a care- 
ful comparision of samples of different quality e.g. 
using a scaling analysis as shown in Fig. 5. Cer- 
tainly, more work is necessary to identify the origin 
of the non-Fcrmi-liquid behavior in systems like 
CeCug-xAua; and to clarify the role of strong dis- 
order. 
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